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Abstract 

The classical integral cross sections of large superfluid ^He^j droplets and 
the number of atoms in the droplets (A^ = 10'^ — 10^) have been measured in 
molecular beam scattering experiments. These measurements are found to be 
in good agreement with the cross sections predicted from density functional 
calculations of the radial density distributions with a 10-90% surface thickness 
of 5.7A. By using a simple model for the density profile of the droplets a 
thickness of about 6-8A is extracted directly from the data. 
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I. INTRODUCTION 



The study of the properties of ^He clusters is currently an active area of theoretical 
and experimental research |l|J^. The interest is based on the fact that '^He clusters provide 
an easily accessible example of a finite size quantum system of strongly correlated particles. 
One of the primary aims of He cluster research is to search for manifestations of superfluidity 
in mesoscopic systems. Recent spectroscopic studies of glyoxal molecules embedded in "^He 
droplets with about 5000 atoms produced in free jet expansions provide the first evidence 
that these droplets are indeed superfiuid The internal temperature of these droplets has 
also been measured spectroscopically to be about 0.4 K in agreement with theory 

Even though there has been notable progress in theoretical work on He clusters during 
the last years the agreement between theory and experiments is not yet entirely satisfactory 
The implementation of ah initio calculations still remains difficult for inhomogeneous 
systems with more than a few hundred atoms, while phenomenological theories, which work 
well for describing macroscopic properties, are often not adequate for systems of finite size. 
A major progress to fill this gap has been made in the development and use of density 
functional (DF) theory p. 

The understanding of the density profiles in the surface region of a quantum fiuid has 
long been considered a very important basic problem [^j. For this reason significant effort 
has gone into the experimental determination of the thickness of the bulk liquid-vapor ^He 
interface 0. At least three different experimental methods have been used to measure the 
surface thickness of bulk liquid He. The first involved atomic scattering experiments p| 
yielding an interfacial width (10-90%) of about 4A. A surface thickness of about 9A was 
measured with an ellipsometric method 0. Finally, x-ray measurements give a surface 
thickness of the order of SA, providing also information on the shape of the interface [llO 



In droplets the surface atoms represent a significant fraction of the total number of atoms. 
For example, in a droplet with 5000 atoms about one half of the atoms are located in the 
surface region. Thus the physical properties of "^He droplets are strongly infiuenced by the 
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surface. Moreover, the surface region is particularly interesting because recent calculations 
predict that the Bose-condensed fraction, which in bulk "^He is about 10%, approaches 100% 
in regions where the density is still within 10% of the bulk density [pJ]] . 

In this paper the liquid-vapor interface of large ^He droplets (A^ = 10^ — 10^) is inves- 
tigated by a new method. The experiment consists in measuring both the integral cross 
section of the droplets and also their number size distributions. From the integral cross 
section and the average number of atoms an effective droplet volume as well as an average 
density is determined. The results are compared with theoretical predictions. Although a 
large number of calculations have been performed for the density profiles of small ^He clus- 
ters (A^ < 10^) |ll|JT^-|T8[ no calculations for the large droplets investigated here have been 



reported so far. In the present work DF theory is used to calculate the density profiles of 
large ^He droplets. The density functional used was recently introduced by the Orsay-Trento 
collaboration [Q. This functional provides density profiles and energies [0 very close to 



the results of ab initio Monte Carlo calculations p!3|-[T5|] in the case of small clusters and 
has the advantage that it can easily be applied to droplets with several thousands of atoms. 
These density profiles provide for a direct comparison with the experimental data and good 
agreement is found within the experimental errors. 

The droplet beam scattering apparatus is described in the following section. The proce- 
dure used to measure the number of atoms in the droplets and their classical integral cross 
sections are described in Section P A| and |IIB| . In Section P^I C| the experimental results are 
presented and in Section |ITl| they are compared with the DF predictions. A final discussion 



is given in Section IV. 



II. EXPERIMENTAL 

Since the experimental methods and the apparatus used in this work have already been 
described in some detail in previous publications |T9| , |20| the experimental procedure will be 



only briefly described here. The He droplet beam is produced by a free jet expansion of 
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^He from a stagnation source chamber at a high pressure and at low temperature (typically 
Pq = 5 — 100 bar, Tq = 4 — 30 K) through a thin walled (20 ^m) nozzle with 5 ± 1 fim 
diameter. After the droplets have passed the skimmer and several differential pumping stages 
they are scattered by a secondary beam produced in another free jet expansion beam source 
which crosses the droplet beam at an angle of 40°. The scattered droplets are detected with 
a homemade electron impact ionizer optimized for a high ionization efficiency followed by a 
magnetic mass spectrometer. The detector can be rotated around the scattering region in the 
plane of the two beams. To avoid capture collisions with the residual gas the background non- 
helium pressures in the differential pumping stages between source and scattering chamber 
were kept below 10"'' mhar. 



A. The Number of Atoms in the Droplets 

The atom number distribution P{N) of the helium droplets is determined from the 
angular distribution resulting from scattering by the secondary beam. A small portion (5- 
10%) of the droplets is deflected by the momentum transfer imparted by single collisions of 
the secondary beam gas atoms, most of which are captured by the droplets . 

The measurements described in the present study have been carried out at source pres- 
sures of Pq = 40 bar and Pq = 80 bar and nozzle temperatures between Tq = 13 K and 
To = 26 K. For these conditions the droplet velocity distributions are very narrow (^ ^ 2%) 
with well known mean speeds I^TJ. The secondary beam gases Ar and Kr were expanded 
from stagnation temperatures of about 300 K and several hundred millibar of stagnation 
pressure from a 40 fim diameter nozzle. The secondary beam contains no appreciable amount 
of dimers or larger clusters, and has a narrow velocity distribution (^ ^ 20%). 

The angular distribution of the deflected droplet beam is measured by rotating the 
detector in small angular increments (typically 300 firad) around the scattering center. Due 
to the large mass of the droplets the deflection angles are very small and a high angular 
resolution is necessary. This is achieved by coUimating the incident beam with a 50 /xm slit 
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in front of the scattering center and a 25 fim slit in front of the detector. The width of 
the expansion zone and the broadening due to colhsions with the residual background gas 
(mostly helium) leads to a FWHM of the beam profile of about 1.5 mrad which defines the 
effective angular resolution. 

Fig. shows three typical measured angular distributions with and without a flag in front 
of the secondary beam [jl9|. The small difference between the two signals is attributed to 



droplets which were deflected after capturing a secondary beam atom. Assuming complete 
momentum transfer the angle of deflection is directly dependent on the number of atoms 
in droplet. As is discussed in more detail below, the present experiments give further jus- 
tiflcation for the hypothesis of complete momentum transfer. The measured droplet atom 
number distributions can be very well fitted with a log- normal-distribution |lT9|,pO|: 



{\nN-fj,y 

where the mean number of atoms and the width (FWHM) are 



(1) 



N = exp(^fi+^-pj (2) 

and 

AA^i = exp (/X - 5^ + 5V21n2) - exp (/x - 5^ - 5^2 In 2) , (3) 

respectively. The FWHM is comparable with A^. Table |I| lists the results for A^, A A"! , S 
and /i , measured using Ar and Kr as defiecting atoms for a wide range of different ^He 

source pressures and temperatures. The values of A^ and AA^i obtained using Ar or Kr for 
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identical "^He source conditions do not always agree exactly, since they depend sensitively 
on the nozzle-diameter and on the nozzle-skimmer distance which was optimized for each 
experiment and differed slightly. In particular the mean droplet sizes for Pq = 80 bar 
with krypton as secondary gas are somewhat smaller by about 30% since the nozzle used 
in these experiments had a diameter of only about 4 /im as estimated from the total gas 
flow into the nozzle chamber. These deviations however have no appreciable effect on the 



average densities. The latter depend on the number of atoms in the droplets and the integral 
cross sections, both of which are determined in the same experiment under identical source 
conditions (see below). 



B. Integral Cross Sections 

The integral cross sections of the droplets are determined by measuring the attenuation 
of the forward peak in the deflection pattern (see Fig. [l|). The attenuation is related to the 



integral cross section a of the droplets according to Beer's law ||20|| , 

— = exp 

-'0 



(4) 



^drop 

where / and /q are the intensities of the ^He droplet beam with and without attenuation, 
Usee is the number density of the secondary beam gas atoms in the scattering center and Lgg 
is the effective length of the scattering volume. Moreover, fdmp is the speed of the droplets 
and frei the relative collision velocity. F^o takes account of the smearing resulting from the 
velocity distributions of the two nozzle beams and leads to a correction smaller than 
about 1%. 

The product of the density in the scattering center and the effective length of the scatter- 
ing volume {usec ■ L^q) was calibrated to within an error of approximately 5% by measuring 
the attenuation of a nearly monoenergetic helium atomic beam, for which the integral cross 
section can be accurately calculated quantum mechanically from the well established inter- 
action potential p3[. The density was calibrated three times with krypton and two times 
with argon at different collision energies and there was no evidence of systematic errors. The 
values of the absolute integral cross sections of the droplets are also listed in Table |. 

The interpretation of the integral cross sections in Eq. (^ is straightforward only in the 
ideal case of an infinite angular resolution. In this case one gets the quantum mechanical 
integral cross section, which includes the forward peaked diffraction part. In the present case 
of finite angular resolution the measured integral cross section is smaller since a fraction of the 
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forward diffraction part of the differential cross section is not included. In order to estimate 
this correction, the quantum mechanical differential cross section for elastic scattering was 



calculated using a partial wave expansion method |^^. The calculations, which are described 
in detail in Appendix indicate that the contribution of the diffraction part to the elastic 
cross section is negligible with the angular resolution of our apparatus. Thus the measured 
integral cross section cr™p turns out to be very close to the classical cross section (Xciass- The 
two cross sections agree to within 1% for droplets with 1000 atoms and for larger droplets 
the agreement is even better. 

The finite width of the beam-profile also has the effect that the deflections of the largest 
droplets in the atom number distributions are too small to lead to a measurable attenuation 
of the beam. This error however can be estimated since the atom number distributions of 
the droplets and the beam-profiles S{Q) are both well known from the measurements. The 
calculation yields a correction of about 1% for mean atom numbers of 10^ and of about 10% 
for droplets with 10^ atoms. 



C. Mean Droplet Densities 

Measured values of the integral cross section, a, as a function of the mean atom numbers 
are reported in Fig. ^. The values obtained with argon and krypton as secondary gas fall 
on a common curve and so we can conclude that the integral cross section is independent 
of the nature of the secondary beam, within the estimated accuracy. The overall experi- 
mental errors are estimated to be about 5% and result mainly from the uncertainty in the 
determination of the density of the secondary beam atoms in the scattering volume. 

From the measured cross section an effective mean density p, defined as the density of a 
uniform sphere with a sharp step edge (liquid drop model) having the same classical integral 
cross section, is determined 



The values of p, normalized to the bulk helium density pbuik = 0.0218 A~^, are given in 
the last column of Table | and the integral cross sections for spheres of different density 
are shown in Fig. § as dashed lines. Here the average integral cross section a is calculated 
from the measured log-normal atom number distributions P{N) by means of the following 
equation: 

2 

a{N) = J^ P(iV) 7r3 i^j dN . (6) 

As can be seen, the experimental effective droplet density goes from about 0.5 ■ pbuik, for 
droplets with 10'^ atoms, to 0.7 ■ Pbuik for droplets with 10^ atoms. This trend is consistent 
with the fact that a significant part of the atoms are in the outer surface region, where the 
density is less than the bulk density. 



III. COMPARISON WITH THEORY 

For comparison with the above experimental results the integral cross sections were 
calculated using density functional theory. These calculations were carried out with the 
improved density functional recently developed |^ for the accurate calculation of ^He droplet 
properties. The energy of the system is assumed to be a functional of the complex function 
^: 

E = Jdrn[^,^*] , (7) 

vl/(r,t) = $(r,t)expQ5(r,t)) . (8) 

The real function $ is related to the diagonal one-body density by p = "l*^, while the phase 
S fixes the velocity through the relation v = {l/m)VS. In the calculation of the ground 
state, only states with zero velocity are considered, so that the energy is simply a functional 
of the one body density p(r). A natural representation is given by 

E = jdvn[p\ = EM + /^r|^(^v^)' ' (9) 
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where the second term on the r.h.s. is a quantum pressure, corresponding to the kinetic 
energy of a Bose gas of non- uniform density. The quantity Ec[p\ is a correlation energy, which 
incorporates the effects of dynamic correlations resulting from the interactions between the 
individual He atoms. The equilibrium configurations are obtained by minimizing the energy 
with respect to the density. This leads to the Hartree-type equation 



where ?7[p, r] = 6Ec/Sp{r) acts as a mean field, while the chemical potential /i4 is introduced 
in order to ensure the proper normalization of the density to a fixed number of particles A^. 
The density dependence of the correlation energy is parameterized in a phenomenological 
way, by choosing a functional form compatible with basic physical requirements and fixing 
a few parameters in accordance with the known properties of the bulk liquid. The detailed 
form of Ec is given in Appendix |B[ The DF theory is particularly suitable for the calculation 
of the density profiles of relatively large droplets. In fact, the density functional of Ref. 
is accurate, by construction, in the limit of the uniform liquid. It also has been tested in the 
opposite limit of small clusters with 20 - 70 atoms |]T^ for which it provides density profiles 
in close agreement with ab initio Monte Carlo calculations. 

The density profiles shown in the upper part of Fig. ^ are obtained by solving Eq. ([To|) for 
= 10'^ to 10^. In order to estimate the contribution to the effective integral cross section 
(Teff from the outer region of the theoretical profiles, the transmission coefficient T{b) of the 
droplets is calculated as a function of the impact parameter b for argon and krypton using 
Beer's law (see the inset in Fig. ^d). For this the trajectory of relative motion is integrated 
along a straight path z, which is a good approximation in the outer region where the density 
is about 1% of the core density (see below): 



The atom-atom integral cross sections a{Ej.e\) were again calculated using the method given 
in Appendix ^, with the relative velocities of the experiment. Here, the more realistic 




(10) 




(11) 
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Tang-Toennies potentials p3| are used in place of the Lennard-Jones(12,6) potentials. The 
calculated transmission T{b) of the droplets is illustrated in Fig. The differences between 
argon and krypton are less than 1% and can therefore be neglected. The transmission for 
densities larger than about 10% of the central density is almost zero, so that the cross section 
is largely determined by the outer region of the droplet. Since the transmission rises very 
steeply the effective droplet-border is rather sharp and corresponds to a radius (Res) where 
the density has fallen to approximately 1%. The integral cross section of the droplets is then 
calculated with the following equation: 

POO 

cTeff = 27r / {1-T{b))bdb . (12) 
Jo 

In a subsequent step the cross sections are averaged over the distribution in the number of 



atoms P{N). The calculated results for the effective radii Rcs = ^Jo'cs/'^, o"efr, the average 
relative densities p/ pbuik and the 10-90% thickness t are reported in Table The average 
relative densities p/ Pbuik in Fig. ^ (filled diamonds) lie between 0.64 and 0.8 and are somewhat 
larger than the experimental values by about 10-25%. This discrepancy suggests that the 
theoretical surface thicknesses are somewhat smaller than the experimental values. 

A simple model calculation can be used to check the effect of the surface thickness on 
the integral cross section. Since all theoretical calculations predict the core density of these 
large droplets to be equal to the bulk density, the density profile can be rather realistically 
described by a simple analytic function [p5| , p6| : 

pW = ^(l-tanh(2l-:?)) . (13) 

where R indicates the point where the density is reduced to 50% of the central density and 
g is a. parameter controlling the surface thickness. By relating the density profile to the 
integral cross section using the same procedure described in Section JTT| an experimental 
surface thickness is estimated by fitting the experimental cross section for the different 
droplets. The resulting values of t are shown in Fig. ^. We obtain a mean 10-90% thickness 
of 6.4 ± 1.3 A, which is somewhat larger than the results of the DF calculation (Table |l|) 
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but still within the error limits. One has to stress here that the results are dependent on 
the model used for p(r); in particular, it assumes a symmetric density profile. Most of the 
predicted density profiles are however slightly asymmetric, with a steeper slope in the outer 
part of the surface. By using an asymmetric profile, with an asymmetry similar to the 
one of the DF profiles in Fig. |^, the estimated surface thickness is about 1-2A larger than 
the values shown in Fig. ^. This would increase the discrepancy between the experimental 
estimates and the theoretical surface thickness. 



IV. DISCUSSION 

As discussed above it appears as if the experimental surface thickness is somewhat larger 
than the values predicted by the DF theory. This is indicated by the 10-25% differences in 
the average densities in Fig. ^ and by the larger widths obtained using model profiles similar 
to Eq.(|13|). Since these differences are greater than the experimental uncertainties of 5%, 
we consider here possible sources of error. 

The possibility that the measured droplet atom numbers are too small, which would 
shift the mean effective density to smaller values, seems unlikely. This would require a 
momentum transfer in the deflection experiment even larger than the momentum Psoc of 
the impinging secondary atom. For elastic backscattering of the atoms in central collisions 
from a hard sphere the momentum transfer would indeed be 2psec- The momentum transfer 
averaged over all impact parameters for elastic scattering from a hard sphere is, however, 
equal to the momentum of the incoming particle Psec- Thus even in the unlikely case of elastic 
recoil scattering, which in view of the large mass of the argon and krypton atom compared 
to that of the He atom seems very unlikely and is not even expected for the scattering from 
a liquid helium surface P, the same results would still hold. Recoil scattering would also 
be in contradiction to the observation that the secondary gas is captured by the He droplets 
I9| , |20[] . Another possibility is a backward directed vaporization of helium induced by the 



impinging secondary gas atom, as proposed by Gspann According to ordinary fluid 
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dynamics this would also seem implausible, since the impact speeds are greater than the 
velocity of first sound of about 240 m/s and therefore a conical Mach shock should be 
created which degenerates in a spherical sound wave. This sound wave would more probable 
induce vaporization in forward direction which would tend to shift our results to greater 
droplet atom numbers N and would therefore increase the discrepancy with the theoretical 
results. 

The comparison between measured cross sections and DF calculations is affected by at 
least two possible sources of problems. The first concerns the form and extent of the density 
profiles, and the second the overall shape of the droplets. The effect of different density 
profiles was already discussed in the previous section and has a noticeable effect. The 10- 
90% surface thickness of the density profiles from the DF calculations is a t = 5.7 A. This 
value is slightly less than other estimates. For example, previous DF calculations for a planar 
free surface p8[ gave t = 7 A, while from ab initio calculations on small clusters [p!^ -[T8[] 



values between 6 and 7 A are predicted. Some experimental and theoretical results for the 
surface thickness of the bulk liquid and droplets are summarized in Table pTT] . 

The second possible reason why the measured cross sections are larger than the calcu- 
lated ones from DF could be that the helium droplets are not spherical, as assumed in the 
analysis, but have the shape of an oblate ellipsoid. A fiattening of the droplets does not 
seem unfounded, due to the large angular momentum of several 1000 h which some of the 
droplets may have as a result of collisions with the residual gas before arriving in the scat- 
tering region. A simple estimate indicates that in view of the large mass of the droplet, the 
resulting rotation speed is slow compared with the colliding secondary beam atoms. Thus 
the effect of a possible deformation on the experimental results can be estimated from the 
calculated geometrical cross section of a rotational ellipsoid averaged over all orientations. 
In this estimate vibrational distortions, which at the low temperature should only have small 
amplitudes, are neglected. 

The radius of the ellipsoid is assumed to have the length a and the radius of the axis 
of rotational symmetry is given by b. For a> b one gets an oblate and for a <b a prolate 
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ellipsoid. The geometrical cross section cTeii under an angle of view uo can be written as (see 
Appendix y): 



The effective cross section is obtained by averaging numerically over all orientations, 

1 /-T 

oiil = - / du sin(u;) creii(u;) . (15) 
z JO 

The resulting mean cross section is shown as a function of the ratio of the radii h/a in Fig. 
scaled to the cross section of a sphere with the same volume. For a real droplet with lower 
density at the surface the effect will be even greater, since the surface area of the ellipsoid 
is larger. To explain the difference of 5-10% between the DF results and the experiments, 
oblate shaped droplets with a value of h/a ~ 0.6 or prolate droplets with h/a ^ 2.5 would 
be required. Oblate shapes with such small average h/a ratios seem unlikely in view of the 
relatively small fraction (^15%) of the incident beam droplets which undergo collisions with 
the residual gas. Thus we conclude that even in the very unlikely case that the droplets are 
somewhat distorted this has little effect on the measured surface thickness. 



By using a combination of scattering techniques (deflection and attenuation), we have 
shown that it is possible to measure the average densities of large Helium droplets. These 
results have been compared with density functional calculations which have been carried out 
for the large droplets with = 10^ to 10^ Helium atoms studied in the experiment. Overall 
agreement can be considered to be satisfactory but on closer examination the predicted 
effective integral cross sections are too small be about 10-20% and the densities too large 
by about 10-25%, which is larger than the errors which are estimated to be about 5%. An 
attempt to fit the experimental data using a symmetric model density profile yields a value 
for the 10-90% thickness of t = 6.4 ± 1.3A which is larger than the density functional value 




(14) 



V. CONCLUSIONS 
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of 5.7 A and is consistent with the differences in the densities. If an asymmetric profile is 
assumed 1-2A larger thicknesses are obtained which would increase the discrepancy. The 
assumptions made in the analysis were critically examined and an explanation for the small 
discrepancy cannot be provided. As discussed in connection with Table |T| the situation is 
similar to the surface of the bulk where experiments also yield thicknesses which are larger 
than predicted by most theories. 

The experimental method has also been recently applied to ^He droplets ||2^ and there 



the average densities relative to the bulk are even smaller than found here for He, a trend 



which is consistent with recent Thomas- Fermi theoretical calculations ||30|| . 

In the future it is conceivable that the present experimental method can be further im- 
proved to provide a sensitive quantitative probe of the outer surface region of large droplets. 
Here it is interesting to speculate what effect the large condensate fraction in this outer 
region |TlI] would have on the cross section. One can also explore the existence of a Landau 



velocity below which the interaction should disappear, as it happens in the bulk liquid at 
57 m/s. Experiments under such conditions are now possible |3l|] and are envisaged in the 
future. 
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APPENDIX A: INTEGRAL CROSS SECTION 

The quantum mechanical differential cross section for elastic scattering from the sec- 



ondary beam gas atoms was calculated using a partial wave expansion method in order 



to estimate the effect of the forward diffraction contribution on our results. The following 
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spherical model potential obtained by integrating the Lennard- Jones 12-6 potential between 
the scattering atom and the He atoms of the droplet was assumed p2| : 



VN{r) 



(r2 - R^y 



r > R . 



(Al) 



The integration is simplified by assuming the He atoms to be homogeneously distributed 
within a sphere of radius R. Here the effective radius R corresponds to the radius Rn = 
tqN^^^ of a droplet with atoms reduced by the effective radius tq of one atom in the same 
droplet, i.e., R = Rn — tq = (A^^^ — l)ro. Estimates of the effective radius ro of one helium 
atom in a droplet can be extracted from quantum many-body calculations. For instance, 
Pandharipande et al. reported the value ro{N) = 2.24 + O.SSA^"^ + 2.59A^-i. The 
parameters e and n are such that, in the limiting case A^ = 1, the potential V/v reduces to 
the two-body Lennard- Jones He-Kr or He-Ar potential: 



Vi{r) =Ae {K/ry'-{K/r) 







( Rm\^ 












K r ) _ 



(A2) 



where Rm = We use e = 2.67 meV and Rm = 3.70 A for He-Kr, and e = 2.59 meV 

and Rm = 3.40 A for He-Ar [Q. Although the inner repulsive potential is not realistic for 
central collisions with the strongly absorbing liquid core of the He droplet, the potential 
Eq. (|A1|) should be a good approximation for describing the glancing collisions, which are 
the relevant ones in determining the elastic contribution to the integral cross section. 

The integral cross section of real He droplets consists of the elastic cross section cTg 
resulting from diffraction in large impact collisions and the absorption cross section aa in 
more central collisions. In the calculations with the hard core potential Eq. (|Al|) the angularly 
distributed isotropic part of the differential cross section takes account of the effect of the 
absorption cross section a^. In our case, where the de Broglie wavelength of relative motion 
A is much smaller than the droplet radius, all particles with impact parameter smaller than 
approximately Rn are expected to be absorbed. The contribution of particles with angular 
momentum / to the cross section is given by ct; = (2/ + 1)ttX'^ iS^. Then the integral cross 
section is given approximately by: 
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Rn/^ oo 

a ^ E (2^ + 1)^^' + E ^.M ^ + ^e. (A3) 

Z=0 1=Rn/X 

The part CTg was calculated quantum mechanically exactly for the potential Eq. (|A1|) . Since 
the calculated elastic scattering is found to be sharply peaked in the forward direction its 
contribution to the integral cross section can be neglected for center of mass angles larger 
than about 5°. The laboratory scattering angle B is in the case of mdrop ^ ^sec related to 
^ by 



6(^9) = ^''''^ <^ Vrci sin 

^drop l^^sec ~r '^drop ) ( 



arcsm + iy 

V ^rcl 7 



-t;secSin/3> , (A4) 



where f rei is the relative speed of the two beams and f3 is the angle between them. Since the 
mass mdrop of the primary beam droplets is much larger than the mass msec of the secondary 
gas atoms the angle G turns out to be very small. For example a center of mass angle = 5° 
corresponds to laboratory angles G of about 1 mrad for droplets with 1000 atoms, which is 
comparable with the angular resolution of our apparatus. 

The theoretical differential cross sections were convoluted with the angular distributions 
5'(G) of the incident beam (beam profiles). These latter angular distributions were measured 
with the same method as described in Section [11 A| with the secondary beam blocked with a 



beam flag located in the scattering chamber. In this arrangement the residual gas pressure 
and its effect on the beam profile is the same as in the scattering experiment. The predicted 
effective integral cross section is then given by: 

where the "i^app represents the effective geometrical center of mass resolution given by the slit 
in front of the scattering center and in front of the detector and s('i9) is the measured angular 
profile of the unscattered beam transformed into the center of mass system. The quantity 
{da /duj) th is the calculated elastic differential cross section. Fortunately the resulting value 
of (7™p turns out to be very close to the classical cross section cTdass = ^ri?^; the two 
cross sections agree to within 1% for droplets with 1000 atoms and for larger droplets the 
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agreement is even better. Thus the corrections discussed above justify the assumption that 
the classical cross section is measured. 



APPENDIX B: DENSITY FUNCTIONAL 

The explicit expression of the correlation energy Ec, entering the density functional 
is given by 

Ec[p]= jdr{\jdr' p{vm\v-r'\)p{v') + |p(r)(p^(r))2 + |p(r)(p^(r))3 

The two-body interaction Vi is the Lennard- Jones interatomic potential, with the standard 
parameters a = 2.556 A and e = 10.22 K, screened at short distance (V/ = for r < h, 
with h = 2.1903A). The two terms with the parameters C2 = —2.411857 x 10^ K A^ and 
C3 = 1.858496 X 10^ K A^ account phenomenologically for short range correlations between 
atoms. The weighted density Ph(r) is the average of p over a sphere of radius h centered in 
r. The last term, depending on the gradient of the density in different points, is introduced 
to reproduce the static response function in the roton region. The function F is a simple 
Gaussian, F{r) = n-^^^e^^ exp{-r^ /f) with £ = 1 A, while = 54.31 A^ and pos = 0.04 

A-3. 



APPENDIX C: MEAN GEOMETRICAL CROSS SECTION OF A ROTATIONAL 

ELLIPSOID 

A rotational ellipsoid is assumed with the axis of rotational symmetry b parallel to the 
z-axis and the other two radii a: 

^ + ^ = 1. (CD 

This can be written in spherical coordinates: 
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Since the visible cross section (Teu of the elhpsoid from a viewing point with a polar-angle lo 
does not depend on the azimuth-angle the problem can be reduced to the x-z-plane. The 
visible cross section Oew is given by the projection of the cutting plane through the ellipsoid 
which is defined by the tangents with angle oj. Since every cutting plane of an elhpsoid is 
an ellipse, the area of the cutting plane Os can written as: 

(Jsiuj) = 7iars{u!) . (C3) 

The projection on a plane perpendicular to the viewing direction is given by: 

(Teii(w) = (Ts{uj) C0S(C^ - Qs) ■ (C4) 

To find the vector Tg = ) = lying in the cutting plane the equation of a tangent 
at the ellipse has to be calculated 

1 = ^ + ^ 

which yields 

6^ b'^x 



Zg d Zg 



X . (C6) 



The tangent must be parallel to the viewing direction u 



With equation ( |ClD follows 



and for z^ and r. 



cot(..) = ^ = . (C7) 

ax a^z^ 



2 ^ a'cot^(^) , . 

^ b^ + a'^cot^u)' ^ ' 



2 

^ b^ + a? cot^M ' 
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r=x, + z 



2 2 6^ + cot^(a;) 



62 + cot^' 



From Eqs. ([Cq ) and (|C9D 6^ can be calculated: 



tan(6s) = — 



Xs a"^ cot{iu) 

The visible cross section of an ellipsoid follows from Eqs. ( p3|) and 



(TelllCo'l = TT a 



¥ + a^cot^{iu) 



cos 



uj — arctan 



coiiuo] 



The mean cross section is derived be integration over uj 

1 

= o / '^^ sin(t^) (Teii(u;) . 
2 JO 
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TABLES 

TABLE L Experimental results as function of source temperature (Tq), source pressure (Pq) 
and secondary gas. The mean number of atoms N and the half-width AA^i are the results of fitting 

2 

the measured mass distributions from the deflection experiment with a log-normal distribution 
(parameters 5 and fi, see Eq. (|l])). The mean classical integral cross section a is obtained by 
attenuation of the droplet beam with the secondary beam. From this data the mean density of the 



droplets p 


as a fraction 


of the known bulk density (pbulk 


= 0.0218 I 


1-3) is 


obtained directly. 


To [K] 


Po [bar] 


sec. gas 


A T 


2 


6 




— r X 2i 


P/Pbulk 


24.0 


40 


Kr 


703 


667 


0.426 


6.46 


2266 


0.40 


22.0 


40 


Kr 


1700 


1632 


0.407 


7.36 


3138 


0.59 


20.0 


40 


Kr 


2617 


2373 


0.528 


7.73 


4519 


0.53 


18.0 


40 


Kr 


4700 


4158 


0.573 


8.29 


6259 


0.58 


17.0 


40 


Kr 


6130 


5331 


0.603 


8.54 


7108 


0.62 


16.0 


40 


Kr 


7741 


6484 


0.662 


8.74 


7661 


0.70 


15.0 


40 


Kr 


8900 


7719 


0.607 


8.91 


8540 


0.69 


13.5 


40 


Kr 


13000 


11240 


0.612 


9.29 


9538 


0.85 


26.0 ^ 


80 


Kr 


1460 


1298 


0.565 


7.13 


3106 


0.51 


24.0 


80 


Kr 


2700 


2524 


0.468 


7.79 


4111 


0.62 


18.0 


80 


Kr 


5260 


4374 


0.673 


8.34 


6270 


0.65 


26.0 


80 


Ar 


2114 


1967 


0.478 


7.54 


3431 


0.64 


24.0 


80 


Ar 


3103 


2835 


0.514 


7.91 


4443 


0.64 


20.0 


80 


Ar 


6458 


5916 


0.509 


8.64 


6615 


0.70 


18.0 


80 


Ar 


9487 


8530 


0.544 


9.01 


9025 


0.67 



^ The smaller cluster number sizes found with Kr instead of Ar as scattering gas at Pq = 80 bar is 
a consequence of a smaller nozzle with an estimated diameter of about 4/im. 
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TABLE II. The effective radii i?efr = yfyesT^ ^md cross sections a^s obtained using Eq. (12) 
and DF calculations of density profiles. The effective density p as a fraction of the known bulk 
density (/Obuik = 0.0218 A""^) and the 10-90% surface thickness t are listed. 



N 


n r S 1 

-Rcff [A] 


r I 9i 
C^eff[A^] 


P/Pbulk 


t [A] 


1000 


25.75 


2083 


0.64 


5.6 


2000 


31.59 


3135 


0.69 


5.6 


3000 


35.69 


4000 


0.72 


5.7 


4000 


38.94 


4763 


0.74 


5.7 


5000 


41.69 


5459 


0.76 


5 7 


6000 


44.09 


6106 


0.77 


5.7 


7000 


46.22 


6713 


0.78 


5.7 


8000 


48.18 


7292 


0.78 


5.7 


9000 


49.97 


7843 


0.79 


5.7 


10000 


51.63 


8373 


0.80 


5.7 
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TABLE III. The 10-90% surface thickness of ^Hcn droplets and bulk liquid helium from 
different published experimental and theoretical works, where N is the number of atoms in the 
droplets and T the temperature (the values for N = oo refer to the planar free surface). 





AT 

i V 


method 


1 [KJ 


t[AJ 


Osborne (1989) g 


OO 


ellipsometry 


1.8 


9.4 


Lurio et al. (1992) |10| 


OO 


X-Ray 


1.13 


9.2 






extrapol. 





7.0 


Pandharipande et al. (1983) [16 


< ( Zo 


\jr iViU , V iViU 


U 


0.0— (.Z 


btrmgari, Iremer (1987) |1^ 


< 728 


DF 





8.8-9.2 




CXD 


DF 





7 


Sindzingre, Klein, Ceperley (1989) H 


64-128 


PIMC 


0.5-2 


~ 6 


Guirao, Centelles, Barranco et al. (1992) ||2^ 


CXD 


DF 


0-4 


6.5 (at 0.4 K) 


Chin, Krotscheck (1995) [|l| 


20-112 


DMC 





~ 6 




20-1000 


HNC 





~ 6 


Barnett, Whaley (1995) |l| 


< 112 


DMC 





~ 6 
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FIGURES 

FIG. 1. Three typical measured angular distributions for a source pressure of Pq = 80 bar and 
source temperatures of Tq = 25 K (a), 20 K (b) and 17 K (c). Krypton was used as secondary 
beam gas. The signals with (filled circles) and without (open circles) a flag in front of the secondary 
beam are show on a logarithmic scale. The weighted differences of the two signals (diamonds) with 
the standard deviations are shown on a linear scale. The mean number of atoms in the droplets 
are: a) = 2602, b) = 6174, c) N = 9834. The integral cross section is determined from the 
attenuation of the forward peak, ie. mrad. 

FIG. 2. The classical integral cross sections averaged over the measured number size distribu- 
tions as functions of mean number of atoms N . The empty symbols show the experimental results. 
The solid line with filled diamonds is the result of the DF calculations. For comparison, the mean 
classical cross sections of spherical droplets with constant density are indicated as dashed lines 
for different values of the density (yCbuik = 0.0218 A^^). The different empty symbols indicate 
the different experimental parameters: triangle: Pq = 40 bar, Sec.-gas=Kr; square: Pq = 80 bar, 
Sec.-gas=Kr; circle: Pq = 80 bar, Sec-gas = Ar. 

FIG. 3. The density distributions a) calculated with a density functional method for droplets 
with between A^ = 10^ and A^ = 10^ atoms with steps of 10'^. The corresponding transmission for 
a beam of krypton or argon atoms is shown in b). The effective radius i?efr for a droplet with 10^ 
atoms is also indicated. 

FIG. 4. The experimental 10-90% surface thickness of He droplets as function of the mean 
number of atoms A^ assuming the symmetric density profile (Eq. (13)) shown in the inset. The 



mean value of t is 6.4 it 1.3A. The DF values are shown as a dotted line. The different empty 
symbols indicate the different experimental parameters: triangle: Pq = 40 bar, Sec.-gas=Kr; 
square: Pq = 80 bar, Sec.-gas=Kr; circle: Pq = 80 bar. Sec. -gas = Ar. 
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FIG. 5. The effective integral cross section of an ellipsoid averaged over all possible orientations 
in units of the cross section of a sphere with the same volume. The radius b is the axis of symmetry 
and a indicates the other radius. 
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